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Abstract 

This paper introduces a geometric method for proving ergodicity of degenerate 

noise driven stochastic processes. The driving noise is assumed to be an arbitrary 
Levy process with non-degenerate diffusion component (but that may be apphed 
to a single degree of freedom of the system). The geometric conditions are the 
approximate controllability of the process the fact that there exists a point in the 
phase space where the interior of the image of a point via a secondarily randomized 
version of the driving noise is non void. 

The paper applies the method to prove ergodicity of a sliding disk governed by 
Langevin-type equations (a simple stochastic rigid body system). The paper shows 
that a key feature of this Langevin process is that even though the diffusion and 
drift matrices associated to the momentums axe degenerate, the system is still at 
uniform temperature. 



1 Introduction 

This paper is concerned with proving ergodicity of mechanical systems governed by 
Langevin-type equations driven by Levy processes and with a singular diffusion matrix 
applied on the momentums. 

Such systems arise, for instance, when one models stochastically forced mechanical 
systems composed of rigid bodies. In such systems one would like to introduce a certain 
structure to the noise and observe its effect on the dynamics of the system. For instance, 
one would like to apply stochastic forcing to a single degree of freedom and characterize 
the ergodicity of the system. The stochastic process associated to the dynamics of these 
systems is in general only weak Feller and not strong Feller. 

The paper provides a concrete weak Feller (but not strong Feller) stochastic process 
to illustrate this lack of regularity. The example is a simple mechanical system that is 
randomly forced and torqued and that preserves the Gibbs measure. In this case one 
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would like to determine if this Gibbs measure is the unique, invariant measure of the 
system. 

A new strategy based on the introduction of the asymptotically strong Feller property 
has been introduced in [B]. This paper proposes an alternative method based on two 
conditions: weak irreducibility and closure under second randomization of the stochastic 



forcing (see theorem 3.1). Our strategy is in substance similar to the one proposed by 
Meyn and Tweedie for discrete Markov Chains in chapter 7 of [9]. 

Although the Hormander condition ([11] 38.16) can also be used to obtain local reg- 
ularity properties of the semi-group, hence a local strong Feller condition and ergodic 
properties. The alternative approach proposed here doesn't require smooth vector fields 
or manifolds, it can directly be applied to Levy processes and (this is our main moti- 
vation) it allows for an explicit geometric understanding of the mechanisms supporting 
ergodicity. 

For related previous work we refer to [8], [7], [6], [3], [2], [5] and [1] . 

2 General set up. 

Let (^t)tgiR+ be a Markov stochastic process on a (separable) manifold M with model 
space M". 

Let {uJt)o<t be p-dimensional Levy process, i.e. a stochastic process on M*' that has 
has independent increments, is stationary, is stochastically continuous and such that 
(almost surely) trajectories are continuous from the left and with limits from the right. 

We assume that there exists a family deterministic mappings (indexed by < t) 
Ft: M X ([0, t] -^RP) ^ M such that 

\^Xt = Ft-^s{Xs, {uJs' - ^^s)s<s'<t) (2-1) 

Recall that the first three condition defining a Levy process mean that {uJt — ujs)t>s is 
independent of (ws')o<s'<s)) the law oiuJt — uJs depends only on t — s and lim^^o IF'[|'-i^s+i ~ 
wtl > e] = 0. 

Recall also |1H [T2] that since is a Levy process, there exists a 7 G M, a constant 
p X p matrix a, a standard p-dimensional Brownian Motion {Bt)t>o and (At)t>o an 
independent Poisson process of jumps with intensity of measure dt x v{dx) on dt x W 
(such that /]gpmin(l, \zP\)v{dz) < 00) such that 

cut = 7t + (7Bt + Ct + Mt (2.2) 

Where Ct = J2s<t ^s^\^s\>i ^ compound Poisson point process (of jumps of norm 
larger than one) and 



Mt = lim A,le<|As|<i-M ziy{dz)) (2.3) 

•^iO V JzeM.P:e<\z\<l ' 

is a martingale (of small jumps compensated by a linear drift). Recall also that any 



process that can be represented as (2.2) is a p-dimensional Levy process, in particular a 
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p- dimensional Brownian Motion is a Levy process. In this paper, the only assumption 
on the stochastic forcing uj will be the following one: 

Condition 2.1. a is non degenerate (has a non null determinant). 

We will then prove the ergodicity of Xt based on the following geometric conditions 
on F. 

Condition 2.2. Xt is approximately controllable, i.e., for all A,BeM and e > there 
exists t>0 andcpe C°([0,t],RP) so that Ft{A, {(ps - (l)o)o<s<t) eB{B,e). 

This condition is illustrated in Fig. [T] 




A 

Figure 1: Approximate Controllability Condition. The condition states that given A,B £ M 
and e > 0, there exists t > and G C"([0, t], R^) so that Ft{A, {<t)s - <?!>o)o<s<t) G B{B,e). 

Condition 2.3. For all < the mapping (x, (f)) ^ Ft[x, {(j)s — 4'o)o<s<t) is continuous 
with respect to the norm — y|| + 1|</> — ^/^|| where \\cj} — 7p\\ := supo<s<t I'/'s — 00 — (V's~V'o)|• 
Let (/>, . . . , be n+1 deterministic continuous mappings from [0, t] onto M" equal 
to at time 0. For A G M", write 

n 

i^l)o<s<t) (2.4) 

i=l 

Condition 2.4. There exists xq £ M and t > 0, such that in a neighborhood of (xq, 0, 0) .■ 

• (x, (0)o<s<i, A) G{x, (0)o<s<t, A) is differentiable in A. 

• VaG is invertible and uniformly bounded. 

• (VxG)^^ is uniformly bounded. 
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3 Main theorem 



Theorem 3.1. Consider a stochastic process Xt on a manifold M that satisfies con- 
ditions 2.2, 2.3, 2.4' 2.1 and admits an invariant measure fi. Let Pt he the semigroup 
associated to X. Then 

• ^ is ergodic and weakly mixing with respect to Pt. 

• ^ is the unique invariant measure of X. 

Proof. We will need the following two lemmas on the Levy process u. 



Lemma 3.1. Assume that uj satisfies condition 2.1. Let < s < t and (p G C'^([0, t], M^) 

be arbitrary. The laws of [lUs — '^o)o<s<t '^^'^ ('^« — "^o — (0s — 'Po))o<s<t absolutely 
continuous with respect to each other. 



Proof. Lemma 3.1 follows by applying Girsanov's theorem to the diffusive component 



Lemma 3.2. Assume that u satisfies condition 2.1 . Let (j) C^{[0,t],MP). Foralle>0, 
the inequality P[supQ<g<j |0s — — {^s — ^o)\ < e] > holds almost surely. 

Proof. Let e > 0. Let (7, a, z^) be the Levy-Khintchine characteristics of uj. Let rj > 
such that ^ 

z^u{dz) < ^ (3.1) 



^1 



zeRP : 0<\z\<ri 16 

Observe that [12] can be written as 

iot = 7"* + aBt + + Mj> (3.2) 

where 

7^ = 7- / 

:»7<|2|<1 

and 



/ zu{dz), (3.3) 

J zeRp ■.n<\z\<l 



C,'' = ^A,1|A.|>, (3.4) 

s<t 

is a compound Poisson point process (of jumps of norm larger than one) and 

= lim (a,1,<|a.|<^ - t [ ziyidz)) (3.5) 

^i" ^ J zmP : €<\z\<ri ' 

is a martingale (of small jumps compensated by a linear drift). Observe that with 
strictly positive probability exp(— tz^(|2;| > r/)), CI is uniformly equal to over [0,t]. 
Furthermore by the Martingale maximal inequality 

E[sup{(M]')2 : < s < t}] < m(yMlf\ (3.6) 
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and using ([12]) 

E[{M^f] = I zMdz) (3.7) 



'zeRP : 0<\z\<r] 

and Chebyshev's inequality we obtain that 



o<s<t 2 2 

hence 

P[ sup \Ms\ < _ ^ „ 

0<s<t ^ ^ 



sup \M^\>U<\ (3.8) 



sup |Ms| < ^] > 1 - ^. (3.9) 



We conclude the proof of lemma 3.2 by applying Schilder's theorem to Bt and using the 
fact that a is not degenerate. □ 

Let us now prove that /x is ergodic. Let A E B{M) be an invariant set of positive 
yu-measure, i.e., 

PtlA = for every t > 0, fi — a.s. (3.10) 

and fi{A) > 0. We will prove that fJ-{A) = 1. Assume < /u(j4) < 1. Then A"^, which is 
also an invariant set, has strictly positive measure, i.e. fi{A^) > 0. Now let us prove the 
following lemma 

Lemma 3.3. If < fi{A) < 1 then 

• For all y e M and e> 0, fi{A D B{y, e)) > 0. 

• For allyeM and e > 0, fi{A'' n B{y, e)) > 0. 

Proof. We will restrict the proof to A. Since /i(A) > there exists xq > such that 
for all e > 0, fj,{A n B{xo,e)) > (otherwise one would get IJ-{A) = by covering the 
separable manifold M with a countable number of balls such that fi(A n B{x, ex)) = 0). 
Assume that there exists yo £ M and e > such that fi{A n B{yo,e)) = 0. Since 
Xt is weakly controllable (condition 2.2) there exists t > and (j) G C^{[0,t\,W) so 
that Ft{xo, {(ps — <Ao)o<s<t) G S{yo, |). From the continuity condition 
Schilder type lemma 3.2 imply that there exists e' > such that 



2.3 



on F and the 



forah xeB{xo,e'), F[Ft{x, {<Ps - <Po)o<s<t) 13{yo, e)] > 0. (3.11) 

Write Pt the semi-group associated with Xf. Equation ( |3.11 ) leads to a contradiction 
with the fact that 

/ Ptix,A)fi{dx) = fi{A). (3.12) 
J A 

since iJ.{An B{xo,e')) > and for all x G B{xo,e'), Pt{x,A) < 1. 



□ 
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From condition 2.4 there exists xq G M and t, e,a,6,K > and such that for 
X £ B{xo,e), ||(/>||L°°(o,t) < and A G {—6,6)"', (</>)o<s<t, A)_is differentiable in A, 



|VaG| < K and (VG) \ < K. It follows from the condition 2.4 and the continuity 



condition 2.3 that e' € (0,e) can be chosen small enough so that there exists z G M, 
< a' < a, < ez such that for ||(;/'||Loo(o,t) < a' we have for all a,b £ B{xo, e'), 

B{z, ez) C G{a, {4>)o<s<t, {-6, <5)") n G{b, {(t>)o<s<t, {-6, <5)"). (3.13) 



Equation (3.13) is illustrated in Fig. |2j 




Figure 2: Closure Under Second Randomization Condition Illustrated. This condition states 
tliat under a second randomization of tlie noise via A, tiie interior of tlie intersection of tlie range of 
G{A, A) (image of (-5, 5)" by A ^ G(A, .)) and G(B, A) is not void. 



Let T > t. From the previous lemma there exists a G B{xq, €')r\A and b £ i3(xo, e') n 
A"" such that Pr(o, A) = 1 and Pr(&, A"") = 1. Set (X^) to be the process Xt started 
from the point a £ M {b £ M) and set Pa to be the measure of probability associated 
to XI . We obtain from the Markov property that 

¥.[PT-t{Xf,A)\=l and E[PT-t(X',A)] =0. (3.14) 

Write 

n 

■.= Ft{a,{oJs-u:o + Y,hv\)o<s<t) (3.15) 

i=l 



The Girsanov type lemma 3.1 implies that the laws of X" and AT"' are absolutely 
continuous with respect to each other. Hence for all A £ [—6,6)", 

E[PT-t(Xr'\A)] =1 and E[PT-t(x''\^)] =0. (3.16) 

Which leads to 

^-2n I W.[PT-t{X^'^,A)]d\ = l and r^" / ¥.[PT-t{X^t^,A)]d\ = {). 

(3.17) 
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Let 0/ be the event 00(9,4) < oi . Observe that from the Schilder type lemma 



3.2 



the 



measure of probability of ilj is strictly positive. It follows from (|3.17|) and (|3.13p that 



]E[ln,V,; 



PT-t{Xl^^.A)\ dX>0 



(3.18) 



Using the change of variable y = X"'^ we obtain from (3.13) that 

dy 1 



E[lf7, 



a, A I 



(y)- 



> 



Hence 
IE[lf7, 



|VAXM|o(xM)-i(y) dy 

Biz,e^) VAX'^'^|o(X-.A)-l(y)|V;,XM|o(XM)-l(y)- 



(3.19) 



> (3.20) 



We deduce from equation (3.13) and the fact that |v^x°-^|o(X''.^)-ifa) bounded from 
below by K~'^ that 



B(.,..)^'^"*^^'^VAXM|o(XM)-l(y)- 



> 0. 



However a similar computation leads from (3.17) and (3.13) to 



0. 



(3.21) 



(3.22) 



Hence a contradiction. Thus ^ must be ergodic. Let us now prove that /i is the unique 
invariant measure. Assume that fi' ^ fj, is also invariant with respect to the semigroup 
Pt. By the argument presented above /u' is ergodic and it follows from Proposition 3.2.5 
of [To] that fj, and fi' are singular and it is easy to check from the argument presented 
above that this can't be the case (the proof is similar to the one given in theorem 4.2.1 
of |10]). Hence ^ is the unique invariant distribution. The proof of the fact that fi is 
weakly mixing follows from theorem 3.4.1 of [lOj and is similar to the one given at page 
44 of [ID] (theorem 4.2.1). □ 



4 Sliding Disk at Uniform Temperature. 

Consider a disk on a surface as shown in Fig. [3] [Ij . The disk is free to slide and rotate. We 
assume that one rescales position its radius and time by some characteristic frequency 
of rotation or other time-scale. The dimensionless Lagrangian is given by 

L{x,v,d,u;) = ^v^ + ^uj^-U{x) (4.1) 

where v stands for the velocity of the center of mass, lo the angular velocity of the disk 
and fj is a strictly positive dimensionless constant given by a := J /{mr'^) (where r is the 
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Figure 3: Sliding Disk. Consider a sliding disk of radius r that is free to translate and rotate on a 
surface. We assume the disk is in sliding frictional contact with the surface. The configuration space of 
the system is SE(2), but with the surface constraint the configuration space is just R x SO(2). 




Figure 4: Ballistic Pendulum. If the dimensionless potential is U — cos(x), then the sliding disk is 
simply a pendulum in which the bob in the pendulum is replaced by a disk and the pendulum is placed 
within a cylinder as shown. 
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radius of the disk, m is its mass and J its moment of inertia). U : M — > M is an arbitrary 
periodic potential which is assumed to be smooth, and of period one. 

The contact with the surface is modeled using a sliding friction law. For this purpose 
we introduce a symmetric matrix C defined as, 

1 1/cr' 

Observe that C is degenerate since the frictional force is actually applied to only a single 
degree of freedom, and hence, one of its eigenvalues is zero. In addition to friction a 
white noise is applied to the same degree of freedom to which friction is applied. The 
governing stochastic differential equations are 



dx 

de 



vdt 
Lodt 



dv 






dt-cC 


V 


dt + aCV2 


dB^ 


duj 











dB^ 



(4.2) 



where C^/^ is the matrix square root of C. The matrix square root is easily computed 
by diagonalizing C and computing square roots of the diagonal entries (eigenvalues of 
C) as shown: 

Write X := {x,9,v,u}). It easy to check that the Gibbs distribution 



is invariant for (4.2), where (3 



2c/a2, Z 



Z 



-dX 



(4.3) 



/ e-P^dX, and E is the energy of the 



mechanical system and is given by 

E : = 

Define 



Y 



-x + a0 
x + 9 



The system (4.2) can be written 



yi(t) = yi(o) + /o*a.t/(^: 



i2(t) = 1^2(0) -/o^a,[/( 



a+l ■ 



ds 
ds 



c-f{Y2{t) -Y2{0)) + aV2Bt 



(4.4) 



(4.5) 



where 7 = {a + l)/a, a = a{a + l)/\/a'^ + 1 and B := {By + B^ 



'2 is a one dimensional 
Brownian Motion. Observe that condition 2.1 is satisfied with to = B, p = 1 and a = (1). 
Observe also that if is a constant then the quantity —v + au is conserved and the 



system (4.2) can't be ergodic. Let us assume that U is not constant, our purpose is to 



prove that the Gibbs distribution fi is ergodic with respect to the stochastic process X. 
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Remark 4.1. Observe that when U is not constant over a non void open subset of M (say 
i))) Y needs to travel a distance that is uniformly (in e) bounded from below by 
a strictly positive amount to get from (Yi,Y2,Yi,Y2) = (0,0,0,0) to the domain Yi > e. 
It follows that in that situation that the process Y and hence X is not strong Feller and 
theorems requiring this property can't be applied. 

Remark 4.2. Observe also that the condition d'^U 7^ in a neighborhood of xq doesn't 
guarantee that Y is strongly Feller in that neighborhood. For instance observe that 
d'^U{xQ) 7^ and dU{xo) > imply that the drift on Yi is uniformly bounded by a strictly 
positive constant on a neighborhood of (0, ^^xq) it follows that lPe(2/i, y2)[^i < 0] is 
discontinuous in the neighborhood of (0, ^^xq) (e) close to the line yi = 0. 

We believe that the system Y is asymptotically strong Feller so one could in prin- 
ciple obtain the ergodicity of /x by controlling the semi-group associated to Y as it is 
suggested in [6j. We propose an alternative method based on the controllability of the 



ODE associated to Y and theorem 3.1, We believe that it is much simpler to control 
the geometric properties of the ODE associated to X rather than the gradient of its 
semigroup. 

One can also check that the generator of Y satisfies a local Hormander condition 
( [TT] 38.16) at a point xq such that d'^U{xo) 7^ so an alternative method to prove 
ergodicity would be to use that condition to obtain a local regularity of the semi group 
associated to U. Here we propose an alternative method which doesn't require U to be 
smooth and which can be applied with Levy processes. 

Theorem 4.1. Assume that U is not constant. Then the Gibbs measure fj, is ergodic 



and strongly mixing with respect to the stochastic process X (4-2). Furthermore, it is 
the unique invariant distribution of X. 



First let us prove that codition 2.2 is satisfied by X. 



Lemma 4.1. Assume U is not constant. Then Y is approximately controllable. 

Proof. Since U is not constant, there exists ti > such that for tj > ti there exists 
a smooth path Y such that Yi(0) = —xi + aOi, 12(0) = xi + 9i, Yi{0) = —vi + ami, 
1^2(0) = vi + LOi, Yi{ti) = -X2 + 0-02, Yi{ti) = -V2 + 0-UJ2 and 

3m('-^] (4.6) 



Take t2 := ti + ^^^^l^'^}, a u ma and interpolate smoothly Y2 between 12(^1) (ob- 



tained from the control problem (4.6)) and 



(Y2{t2)\_(x2+e2\ 



Observe that the extension of Yi to {ti, t2] as a solution of (4.6) satisfies 

Yi{t2)-Yi{ti] 

^(to) - ^(t 



< \ (4.8) 
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Taking (p be the smooth curve defined by i?i>(0) = and 



-dM 



a + 1 



dZ2 _ 1- 
C7— — + av2 
dt 



dt 



completes the proof. 



(4.9) 



□ 



Proof. The proof that X satisfies condition 2.3 is a standard application of Gronwall's 
lemma. Observe that the semi-group associated to X is not strongly irreducible and 
never equivalent to fi because \{—v + au)){t) — {—v + (Tl<j)(0)| < ||92,.[/||Loot. Let us now 
show that condition 12.41 is satisfied. 

Write ^ the stochastic process defined by 



'ei(t) = ei(o) + /o*a,.f/(^)d. 

Ut) = 6(0) + aV2Bt 



(4.10) 



To prove that Y satisfies condition |2.4| it is sufficient to show that ^ satisfies condition 

[M 

Since U is smooth and not constant, there exists a point G [Oil)i e, C > such 
that for X G B{x^ , e), d'^U > C. Let be a point of the phase space such that 

and Ci = C2 = 0. Let < e' < e/100 and a e -B(C, e'). 

Let , . . . , (/74 be 4 continuous mappings from onto 
For A G we write the solution of 



"•C2-Ci _ ^0 

0-+1 ~ 



equal to zero at time zero. 



m 



<T+1 



ds 



di + /o d,JJ( 

<^i + !l^i{s)ds 
"2 + Jli2{s)ds 



(4.11) 



It follows that 



f)!idM^-^ + aC-^-^-^))is)dsda 



cHt) - eAt) = j^ic^is) - e?(.)) ds 

U2 (*) - ^2it) = E-=i A.(a^/2 y,,{s) ds) 



Writing r] the solution of 

m + /d vis) ds = d^,u{x^) ^ f^ic^ - e2){s) ds 

= lo vis) ds 



(4.12) 



(4.13) 



we obtain that up to the first order in A, and at the order in e' and t, (Ci(i) — 
^i{t)^Ci{'t) - ^lii)) can be approximated by {'f]{t),r]{t)). It follows that - Ct can be 
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written as M(X,t)X where M(X,t) is continuous in t and A in the neighborhood of 0. 
Moreover, ipi, . . . ,(p4, can be chosen so that M, and are uniformly bounded in that 
neighborhood. Choosing < 6 <^ 1 and < e' <^ 6t <^ 1 imphes condition 2.4 By 
invoking theorem 3.1 one obtains that the process is ergodic and weakly mixing. 

It follows from theorem 3.4.1 of [TO] that for (p G L'^ifJ-) there exists a set / C [0, +cxd) 
of relative measure 1 such that 



lim E[{p{xt,9t,vt,ujt)] = fi[(p] in L^(/i). 
\t\-*oo, tei 



(4.14) 



Furthermore since t K[ip(xt, 9t,vt,uJt)] is continuous when (p is continuous and bounded 
we deduce that when (p is continuous and bounded then 



lim E[(f{xt,6t,vt,uJt)] = fJ-[f] in L'^ifJ-)- 



(4.15) 



The fact that the process is strongly mixing then follows from corollary 3.4.3 of |10j. □ 



In [T], using thoerem 4.1 we prove that if U is non-constant then the x-displacement 
of the sliding disk is /i a.s. not ballistic (see Proposition 4.1). However, the mean- 
squared displacement with respect to the invariant law is ballistic (see theorem 4.2). 
More precisely, we show that the squared standard deviation of the x-displacement with 
respect to its noise-average grows like t^. This implies that the process exhibits not only 
ballistic transport but also ballistic diffusion. If U is constant then the squared standard 
deviation of the x-displacement is diffusive (grows like t). See below for theoretical 
results and numerical experiments using efficient stochastic variational integrators. 

Proposition 4.1. Provided that U is non-constant, then fi a.s. 

x{t) - x{0) 



lim 



0. 



Proposition 4.2. The squared standard deviation of the xt + 9t-degree of freedom is 
diffusive, i.e., 

^t,[{xt + et-nxt + et]?] 2aV 



t-^oo t c^{a^ + 1) 

Proposition 4.3. Assume that U is non constant, then 

v2" 



(4.16) 



lim sup 



xt + c7et-E[-xt + c7et]y 



t2 



< 4 



and 



lim inf 



{-xt + aet-E[-xt + a9t]y 



t2 



> 



11 + a 



4 p 



(4.17) 



(4.18) 



Theorem 4.2. We have ('Jl]) 
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// U is constant then 



lim 

i— >oo 



2^2 



c2(cj2 + l)(a + l)2 



(4.19) 



// U is non constant then 



and 



lim sup 

4— >oo 



lim inf 



E^[{xt-E[xt])' 
t^ 

t2 



< 



/3(l + a) 



> 



4/?(l + cj) 



(4.20) 



(4.21) 



Classical homogenization techniques can't be applied to obtain theorem 4.2 (since 
the stochastic forcing is degenerate on momentums). We refer to [1^ for a proof of that 
theorem. The ballistic diffusion is caused by long time memory effects created by the 
degeneracy of the noise and the coupling between the two degrees of freedom through U . 
Figure [5] gives an illustration of the mean-squared displacement of the rolling disk versus 
time started from rest. In [Ij we have used that phenomenon to propose a fluctuation 
driven magnetic motor characterized by ballistic diffusion at uniform. A plot of the 
angular displacement of that magnetic motor versus time for a single realization started 
from rest is given in figure [6] 





Usymmetrto ■ 

- - Uasymmetrto ' 






:,;,;,;^''',;„,„,,v:;:::: 









Figure 5: Sliding Disk at Uniform temperature, h = 0.01, a = 5.0, c = 0.1. A log-log plot of the 
mean squared displacement of the ball. It clearly shows that the x-position exhibits anomalous diffusion 
when f/ is symmetric or asymmetric. The disk is started from rest. In the control and flat U cases the 
diffusion is normal. 
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(a) a = 0.001 (b) a = 0.001 (c) a = 0.00075 (d) a = 0.00075 

Figure 6: Angular position of a magnetic motor (uniform temperature). Four different 
realizations of the angular component of the center of mass of a magnetic motor are plotted. The system 
is started from rest. 
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